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Abstract 

m 
o 



For p = 1 (mod 4), we prove the formula (conjectured by R. 
Chapman) for the determinant of the ^i- x 2±1 matrix C = (Cy) 

with = 



1 Introduction 



Let p be a prime and (-) denote the Legendre symbol. Let us set n = 
and consider the following (n + 1) x (n + 1) matrix C: 

Cij = (~~) ' < z, j < n. 

Here and it what follows it is more natural to enumerate rows and columns 
starting from 0. R. Chapman [8] raised the problem of evaluating det C; for 
motivation and related determinants see also [B], [7]. In particular, Chap- 
man conjectured (see also jH Problem 10]) that detC is always 1 when 



*This research was supported in part by the Dynasty Foundation and RFBR (grant 
09-01-00784-a) and the State Financed Fask project no. 6.38.74.2011 at St. Petersburg 
State University. 

2010 Mathematics Subject Classification: Primary 11C20; Secondary 11R29, 15A15, 
15B05. 

Key words and phrases: Legendre symbol, determinant, Cauchy determinant, Dirich- 
let's class number formula. . 



1 



2 



M. Vsemirnov 



p = 3 (mod 4) and had a conjectural expression for det C in terms of the 
fundamental unit and class number of Q(y/p) for p = 1 (mod 4). The se- 
quence det C for primes p = 1 (mod 4) also appears as sequence A179073 
in the On-line Encyclopedia of Integer Sequences pQ . 

Chapman's conjecture for p = 3 (mod 4) was settled affirmatively in 
[12] . The aim of this paper is to apply methods developed in [12] and to 
evaluate det C for p = 1 (mod 4). 

Let (9 be the ring of integers of Q(y / p). Let e be the fundamental unit 
in O and h = h(p) be the class number. 



Theorem 1. Let 

(1.1) a + by/p 



e h , if p = 1 (mod 8), 
e 3h , ifp = 5 (mod 8) 



T/ien det C = —a. 



Our proof is divided into three steps. First, we decompose C into a prod- 
uct of several matrices; see Theorem [2] below. This part resembles a simi- 
lar step in the evaluation of det C for p = 3 (mod 4); see [12] for details. 
Second, we find general expressions for certain parametric Cauchy-type de- 
terminants and reduce our problem to a particular case of that calculation. 
Finally, we relate the obtained determinants to Diri chief's class number for- 
mula for real quadratic fields [U Ch. 5, §4], which involves both class number 
and the fundamental unit. 

Looking at the numerical data, W. Zudilin and J. Sondow conjectured [H 
Entry A179073] that detC is always negative and even. This easily follows 
from our Theorem 1. 

Corollary 1. For p = 1 (mod 4), detC is negative and even. 



2 Matrix decomposition 

Set n = (p—l)/2. Let ( be the primitive p-th root of unity with arg ( = 27r/p. 
We also fix its square root £ 1 ' 2 such that arg^ 1 / 2 = ir/p. 
Let us consider the following three matrices U, V, and D: 

fiV-i-2* _l (i\f-2j-i 

(2 1) 77" = — < i ?' < n 

(2-1) U v C -^+(i)g) ' °SM<», 



V« = ( 2ij , 0<i,j< 



n. 



(2.2) Da = J] 7577^' < z < 

0<fc<n 

(2.3) A, = 0, / ./. 



On R. Chapman's "evil determinant 



3 



In particular, V is a Vandermonde-type matrix and D is diagonal. If we set 

g(x)= J] {x-C 2 % 

0<k<n 

then the diagonal entries of D can be represented in an alternative way as 

1 



(2.4) A 



Finally, let r p (r) be the Gauss sum 

Theorem 2. For any prime p such that p = 1 (mod 4) ; we /iai>e 

(2.5) C = r p (2)C (p ~ 1)/4 ■ VDUDV = VpC (p " 1)/4 • VDUDV. 

Remark 1. For p = 3 (mod 4) we have (see [12]) a similar expression 

C = -r p (2)C" (p+1)/4 • VDUDV, 

where 

(2.6) LT„ = ^ ^fk • 

Both ( 12. ip for p = 1 (mod 4) and (12. 6p for p = 3 (mod 4) can be unified as 

(;)C-- 2l + (f)C- 2 - 



" — Z-l—l 



n n 1 1 ^"\/"— »•— 2A; _i_ 2) — ft 
2ki+2rj _ 1 ! _ Vp/^ 



The unified formula for the decomposition becomes C = (— )r p (2)^ p 
Proo/ 0/ Theorem^ Let 5 = VDUDV. We have, 

n n 

The term corresponding to k = r = in the above sum vanishes. The 
remaining terms can be arranged into three groups depending on whether 
k = 0, or r = 0, or k ^ 0, r ^ 0. More precisely, 

Bij = s + si + s 2 , 
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where 
so 

Sl 
S2 



1 y^fr\ C (2J_1) 



re n 



2fc) ' 

re n 11 /AA/^— r— 2fc _i_ / r\ f—2r~k 
2ki+2rj ' W ^ W_; 



fc=l r=l 



Notice that (- 2k - 2r ^ 1 and ((J) g)) 2 = 1 for 1 < fc, r < n. Applying the 
identity 

;<-©e))(-e)e))=<--" 

which is valid for /c, r such that p f k, p \ r, we have that 



ra n 



S 2 



EE^7 



c 



(r^-g)©)(e)c- 2fc + ©c- 2 ^) 



^ 9'(C 2k )9'(C 2r ) c- 



■2k-2r 



k 

re n 



1 



n '( 



= S3 + s 4 , 
where 

n re 

» 3 = EE ft 

fc=l r=l 

re ra 

= EE 

fc=l r=l 

re re 

- - EE"' 

fc=l r=l 

re ra 

= EE 

fc=l r=l 

By the Lagrange interpolation formula, for any polynomial / of degree 
less than n + 1, we have that 

/(*) _ y i /(C 2r ) 



^2ki+2rj 


C 


-2r- 


-3, _ C - 


-2r-fc 


g'(( 2k )g'(( 2r ) 




c- 


■2fc-2r _ 


1 


£-2ki+2rj 


C 




-c k 




g'(( 2k )g'(( 2r ) 


c 


-2fc 


-C 2r ' 




£2ki+2rj 


c 


-3r- 


-2k _ f- 


-r-2k 


g'{t 2k )g'{t 2r ) 




C" 


2fc-2r _ 


1 


^2ki+2rj 


c 


-3r 


-c _r 




g'(( 2k )g>(( 2r ) 


c 


-2r 


-C 2fc ' 





Therefore, for any x different from the roots of g, 



/(C 2r ) /(*) 1 /(l) 
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If k runs through 1,. . . , n, then —2k (mod p) runs through odd integers 
p — 2, p — 4,. . . , 3, 1. In particular, (~ 2k is not a root of g. To evaluate 
S3 we make summation with respect to r first and use (12. 7p for f(x) = x^ 
substituting x = (~ 2k : 



S3 



^ w 5 (C 2fe ) 



<?(C~ 2fc ) <7'(1) C~ 2fc 
) 1 A A\ C (2i ~ 1)fc 



" ^2fc(i-j)fc-3fc_^-fc) 

fcj w ^'(c 2fc )^(c- 2fe ) 



Si. 



To evaluate S4 we argue in the same way but now we sum with respect to 
k first and substitute x = (~ 2r into (12. 7ft for f(x) = x\ As a result, 

Therefore, 

Sy = h. w ^'(c 2fc )^(c- 2fc ) + E w ^mc- 2 *-) 



fc=i 



fc\ (£2fc(i-j) _|_ ^-2fe(i-j)V^- 3fc _ £-fc\ 



2fc^i 



Now we evaluate the denominator of each term. Recall that n — (p — l)/2, 
so (-l) n+1 = (-1)(p+ 1 )/ 2 = -1 if p = 1 (mod 4). We have 



g\e k )g{C 2k )= WiC-C) \[{C 2k -C) 

0<t<n 0<t<n 
= ( _ ir+ l (C -2, r+ l (c 2 )( + l + ... + n) JJ( C »_ C ») J] ( C »_ C - 



-2t\ 



0<t<n 0<<<n 



<- fc(p+i) c (p2 - i)/4 (c 2fc - 1) n (c 2fc -n 



0<t<p-l 



= -rt (p - 1)/4 (c 2fe -i)((^-i) / )i,=c- 

= -K- 3fc C (p - 1)/4 (C 2fc -i). 

Using this together with the fact that (-) is an even character for p = 1 
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(mod 4), we continue as follows: 

£-(p-l)/4 « / k 



/4 » /lA 

% = ^ W-V^ + c 

p ^ w 



-2k(i-j)\ 



V \\Pj \ P 



-2fc(i-j) 



-2k(i-j) 



P / P VP/ P 

Since r p (2) = (|)\/p for p = 1 (mod 4) (e.g., see [9J Ch. 6]), we conclude 
that 

c^Vp(2)s« = Q)^ 1 ^ = (^) = 

which completes the proof. □ 



3 Determinants of Cauchy-like matrices 

The following identity is due to Cauchy [5]: 

det ( — - J = JJ JJ (u i + w i )- 1 . 

V » ~r j / i,j=0,...,m-l 0<i<j<m-l 0<i,j<m-l 

An alternative form of this identity can be obtained by replacing Vj with 
vj 1 and multiplying each column by v J : 

(3.1) det ( — —) 

\ 1 + U i V jJi, j =0,...,m-l 

n i( u i - u j)( v i - ^ n c^+uiVj)' 1 . 

0<i<j<m-l 0<i,j<m—l 

For further connections of ( 13. ID with representation theory and symmetric 
functions, see [HI Ch. 7]. 

In this section we evaluate determinants of several parametric matrices 
related to the second form of the Cauchy identity. Let it = (u , . . . , w m _i), 
it = (t> , . . . , f m _i). Assume further that 1 + u(Uj ^ 0, i, j — . . . , m — 1. 
Let M m (Tt, it) be the m x m matrix with 
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Theorem 3. We hav£\ det M m (lt, ~&) = 

_ ,m— 1 m— 1 

= o f it it a + «i) + n c 1 - ^ n ( x 



m— 1 



m— 1 



x i=0 

x n i 

0<i<j'<m-l 



3=0 



Ui — Uj 



i=0 j=0 
m—1 m—1 

n - ^) n n + 

0<i<i<m-l i=0 j=0 



- Vh 



u i v j/ 



Proof. Let J be the m x m matrix with all entries equal to 1. Consider 

f(t) = det(tJ + M m (~J,^))- 

Since J has rank 1, there are two invertible matrices H\ and H2 with coef- 
ficients independent of the variable t, such that 



H1JH2 



(\ ... o\ 
... 

\0 ... oy 



Since f(t) = det tfr 1 det(t#i JH 2 + HxM m {lt, ~J)H 2 ) det#2~\ the func- 
tion / is linear with respect to £. In particular, detM m («, v) = /(0) = 
(/(l) + /(-l))/2. On the other hand, 



/(l) = det 1 + 



Ui + 

1 + Ujfj/ j J= o,...m-l 



det 



(1+^)(1 + ^) 

1 + UiVj 



i,j=0,...,m—l 



m—1 m—1 1 x 

II<i + *)II(i + «a)-<tot(T+^;.. o t - 

i=0 i'=0 v 1 3 J i,j=0,...,m-l 



In a similar way 



/(_!) = det -1 + 



Ui + Uj 

1 + UjVi , , 
3 / t,j=0,...m—l 



det 



-(1 - Uj){l - Vj) 

1 + UiVj , ,. . 

J / %,j=D,...,m—l 

m—1 m—1 



m—L m—i / 1 \ 

= Mrn(i-.)n(i-",)'det( TT ^) . 

i=0 j=0 v * 3/ i,J=0,...,m-l 

Combining this with ( 13. ip we complete the proof. 



□ 



When the paper was ready, T. Amdeberhan informed the author that an equivalent 
statement had been proved in [3] by a different method. 
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Now let it — (xi, ... , x m ), V = (Vi, • • • i J/m) an d assume that 1 + x^yj ^ 
0, 1 + Xi 7^ 0, 1 + yj 7^ 0, z, j = 1, . . . , m. Let W m (x, y) be the following 
(m + 1) x (m + 1) matrix: 



/0 



i \ 



Theorem 4. VFe /lave det W m (x, y) = 

1 , m m m m \ 

= - 2 ( ri( i + ^) ri( i + vi) - ri( i - *o ri( i - ) 

M=l 7=1 i=l 1=1 ' 



x yi i x i~ x i) n - f<) n IB 1 + l- 

l<i<j<m l<i<J<?ri i = l J=l 



Proof. Since 



/ 1 



M m+1 ((l,^),(l,^)) 



1 \ 



M m (~x>,f) 



we have that 

(3.2) detW m (^,^) 



/ 1 



detM m+ i((l,-^),(l, V)) - det 



V 



J 



= detM m+ i((l,^),(l, V)) -detM m (^,^). 
By Theorem [3], 

det M m+1 ((l, (1, )) = (1 1} f[(l + x t ) f[(l + Vj ) 



i=i 



xJJ(l-^) Yl ^i-^YliVs- 1 ) n 

j'=l l<i<j<m j=l l<i<j'<m. 

m m mm 

x n^ 1 + n^ 1 + ^ n n^ 1 + + ° 

j=l j=l i=l J=l 



i=l 
m 



j=l l<*<i<™ l<i<j'<m 



m m 

xiniii+xM)- 1 . 

i=l j=l 



On R. Chapman's "evil determinant 



9 



Applying Theorem [3] to the evaluation of det M m (~5", ~$) and using (13. 2p we 
complete the proof. □ 

4 Evaluation of the determinant 

Recall that n — (p — l)/2. Let G be the diagonal matrix with 

(4.1) G 00 = l, G «=(~V, i = l,-..,n. 
Set 

(4.2) W = GC/G. 
By a direct computation, Woo = 0, 

W 0j = W j0 = l, j = l,...,n, 

(ffc + (gc j . . 



In particular, = W n (2?", y), where Xi = yi = (~)0j i — 1, ■ ■ ■ ,n. Now 
we apply Theorem HI The last product in the expression for det W n can be 
transformed in the following way: first we extract terms with i — j and then 
we join together the two factors corresponding to and (j, i) for i ^ j. 
Taking into account that (— l) n = (— l)^ 1 )/ 2 = 1 for p = 1 (mod 4), we 
conclude that 

«, ^.=-t^(n(,.g),)--n(,-g),'- 
',n (©<-©<')■ n ('*6)gH>«* 



The following auxiliary result is an easy consequence of standard meth- 
ods of evaluating Gauss sums. We present its proof for the sake of complete- 
ness. 

Lemma 1. If p \ r, then 

(4.4) ncc*' - r**) = Qvp- 

Proof. By [9l Proposition 6.4.3], 



= flcc^-r^) = (-i) n ncr-^-c- 



^ -p+(2fc-l)N 
fc=l fc=l 



II(c 2i -r 2i ). 
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Let us apply the automorphism 7 of Q(C) induced by 7(C) = C- It follows 
from the standard properties of Gauss sums that j(y/p) = (^)y/p, which 
completes the proof. □ 

Corollary 2. We have 

(4-5) f[(C j/2 - C m ) =(l)^= ("l) n/2 ^ 

(4.6) f[(l + ^) = C in+1)/2 (^). 

Proof. We have C 1/2 = -( {p+1)/2 . Hence, 

n n 

U(c j/2 - c i/2 ) = Il((-c (p+1)/2 ) J - (-c (p+1)/2 n 



i=i i=i 



= (_i)™/ 2 ^((c^ 1 ^ 2 )^ - ((( p+1 )/ 2 )^'). 

(Here we extract —1 from each term that corresponds to an odd j.) The 
first desired identity now follows from Lemma [1] applied to r = (1 — p)/4 and 
from the fact that, for p = 1 (mod 4), (-l) n / 2 = (-1)(p~ 1 )/ 4 = g). To prove 

the second identity we use the equality 1 + ( 2j = C J (C 2j ~ (~~ 2j )i( j ~ C~ J ) -1 
and Lemma [1] applied to r = 1 and r=(p+l)/2. □ 

Lemma 2. We /iai>e 



5 (n + g) c V - n f 1 - g) c j ) 2 ) = (-1)^^, 



where a is defined in (11. ip . 
Proof. Let 

n(i+(^)c 



s = 



n / / .\ \ 2 



Notice that (1 + (j)^') 2 = (C j + (j)) 2 - We have 

* = n {<> + g)) 2 = c"'" +i,/2 rj (c j/2 + g)c- j/2 ) : 
= c (,,+i,/2 n K j/2 -c- i/2 ) 2 n <c j/2 + c- j/2 ) 2 



1 < j < n 1 < j < n 

07p)=-i 07p)=i 



C n(n + l)/2 "Q (C i/ 2 _ r i/ 2)a "Q ( C ^_ C -i/ 2) 



-2 



1 < j < n 1 < j < n 

07p)=-i 07p)=i 



n n 

n (c j - c 3 ) 2 ]> 2 - c-' 72 )- 1 n ^ 2 - c~ j/2 )- 



l<j<n j = l j = l 

(j'/p)=i 
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Since (-) is an even character for p = 1 (mod 4), the number of quadratic 
residues modulo p on the interval [l,n] equals the number of quadratic 
non- residues on the same interval. Therefore, we can continue as follows 



■\2 / \ -2 / o \ 2 

^3 i TT I ■ n J \ TT I ■ ^3 



n n n 



sin 



l<j <n x 1 < j < n 1< j < n 

07p)=-i 07p)=i U/p)=i 



x 



n / .x -1 n 

n(-7j nK" 2 -c- j/2 ) 



Let us consider the third and the fourth products more carefully. We have 
sin(27r/'/p) = sin(7r(j9 — 2j)/p). Moreover, the map j h-> p — 2j gives a 
bijection between the sets 

{j : n/2 < j < n, (J) = 1} and {A; : 1 < fc < n, fc is odd, g) = g)}, 
while the map j h- >■ 2j is the bijection between 

{j : 1 < J < (j) = 1} and {k : 1 < k < n, k is even, (f) = g)}. 

Therefore, 

sin vj n( sin -j = n sin — n ( sin v) • 

, , ' j' = l ^ ' l<k<n l<k<n ^ ' 

(j/p)=l " (fc/p)=(2/p) (fc/p)=-(2/p) 

By Dirichet's class number formula for real quadratic fields (e.g., see [4, 
Ch. 5, §4]), 

n -7 n H)" 

l<j<n l<j<n ^ " ' 

07p)=-i U/p)=i 
Combining with the above equalities, we conclude that 

In a similar way, 

" //-NX 2 

( 4 - 8 ) II ( 1 ~ ( ~ ) ^ I = C n{n+1)/2 e~ {2 ~ {2/p))h Y[(( j/2 ~ C j/2 )- 
j=i ^ ^ 7 ' 3=1 

It is well known that e h is a quadratic unit of norm —1 (which is equiv- 
alent to the fact that the norm of e is — 1 and h is odd); e.g. see pUl Ch. 4, 
§18.4] or [4, Ch. 5, Sec. 4, Ex. 5]. It follows that 

(4.9) e (2-(2/p))A _ £ -(2-{2/p))ft = 2a? 

where a is defined in (11.11) . 

Applying Corollary [2] to the evaluation of YYj=i(C^ 2 ~C~ 2 ) we complete 
the proof. □ 



3=1 v vPy 7 3=1 
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Proof of Theorem [IJ It follows from Theorem [21 Lemma |2] and equations 

( P , f faj) . ( BID that 



n+2 



(4.10) det C = _ C («+i)(^i)/4 ^ Q ^ 

x (det V) 2 (det D) 2 (det G)~ 2 flf^ 2 a, 

where 

A - n (Q- (;V 



l<i<j'<n 



/. - n ^ c 

and a is defined by ([TT]) . Clearly, det G G (Z[C])* by (|Q]i . 

Let us recall some well-known facts about the arithmetic of the ring 
Z[£]. The reader may find further details in Ch. 13]. The ideal (1 — £) is 
prime in Z[£] and p = a(l — C) p ~\ where a is a unit in Z[£]. In particular, 
^ = r p (l) = a x (l - C) n , where a, G (Z[C])*. Finally, (1 - C c )/(1 - C) and 
1 + ( c are in (Z[£])* provided p \ c. 

Notice that 

detv = n (c 2j -C 21 ) 

0<i<j<n 

by the well-known evaluation of the Vandermonde determinant. Using this 
together with the definition of D (see (12. 2p . (12.31) ) and the above observa- 
tions we find that 



( v ^) n+1 (dety) 2 (det J D) 2 G (Z[ 

Since (-) is an even character and n = (p— l)/2, there are n/2 quadratic 
residues and n/2 quadratic non-residues modulo p on the interval [l,n]. 
Thus, 

#{(iJ)-l<iJ<n> g) = (J)} = ^/2, 
#{(i,i):l<i,J<n, (J) =-(J)} = n 2 /2. 

Since the pairs (i, i) are in the first set and the pairs and (j, i) are in 
one and the same set, 

: 1 < i < j < n, g) = (J)} = n(n - 2)/4, 
#{(i,j):l<i<j<n, (i) =-g)} = n 2 /4. 

In addition, 

^ \ (1-C)(Z[CD*, if © = ©, 

i + nr^i e I (z[c]r ' = 

WWi 1 (i-C)(z[CD*, ^ ©^©- 
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Therefore, y/pfifc 2 G (Z[C])*. 

Finally, notice that ( is of odd multiplicative order and therefore any 
power of ( is a square in Z[Q. Using (I4.10p we conclude that det C = —a5 2 , 
where 5 G (Z[£])*. Since det C is an integer and a is a non-zero integer 
or half-integer, we have that 5 2 G Q. Hence, 5 2 G Z*, i.e., 5 2 = ±1. On 
the other hand, Q(C) does not contain primitive fourth roots of 1, since 
[Q(C, ^T) : Q] = [Q(C • <fi) : Q] = 2(p - 1) and [Q(C) : Q] = p - 1. 
Therefore, 5 2 = 1. This completes the proof. □ 

Proof of Corollary QJ Since £ > 1, we have immediately that a > 0. Hence, 
det C < 0. 

Now write e ft = (a + /3 v /p)/2, where a and /3 are integers. Since e h is a 
quadratic unit of norm —1 (see [TOi Ch. 4, §18.4] or [4, Ch. 5, Sec. 4, Ex. 5]), 
we have 

(4.11) a 2 -p/3 2 = -4. 

Consider (14.111) modulo 16. A direct search shows that 

• for p = 1,9 (mod 16), a = 0, 4, 8, or 12 (mod 16), i.e., a = a/2 is 
even; 

• for p = 5,13 (mod 16), a 3 + 3a/3 2 p = (mod 16). In particular, 
raising to the third power, we have that a = (a 3 + 3ct/3 2 p)/8 is 
even. 

□ 
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